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ABSTRACT 

The nature of ultraluminous X-ray sources in nearby galaxies is one of the 
major open questions in modern X-ray astrophysics. One possible explanation 
for these objects is an inhomogeneous, radiation dominated accretion disk around 
a ~ 1OM black hole — the so-called "photon bubble" model. While previous 
studies of this model have focused primarily on its radiation-hydrodynamics as- 
pects, in this paper, we provide an analysis of its X-ray spectral (continuum 
and possible edge and line) characteristics. Compton reflection between high 
and low density regions in the disk may provide the key to distinguishing this 
model from others, such as accretion onto an intermediate mass black hole. We 
couple a Monte Carlo/Fokker-Planck radiation transport code with the XSTAR 
code for reflection to simulate the photon spectra produced in a photon bubble 
model for ULXs. We find that reflection components tend to be very weak and in 
most cases not observable, and make predictions for the shape of the high-energy 
Comptonizing spectra. In many cases the Comptonization dominates the spec- 
tra even down to ~ a few keV. In one simulation, a ~ 9 keV feature was found, 
which may be considered a signature of photon bubbles in ULXs; furthermore, 
we make predictions of high energy power-laws which may be observed by future 
instruments. 



Subject headings: accretion, accretion disks — methods: numerical — radiative 
transfer — X-rays: binaries 
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Introduction 



Approximately 150 off-nuclear Ultraluminous X-ray Sources (ULXs) in nearby galaxies 
have been discovered with luminosities great er than 10 39 erg s -1 , exceeding the Edding- 



ton Luminosity for a ~ lOM^ black hole (e .g., Fabbiano 19881 ; Colbert &: Mushotzky 1999 



Colbert fc Ptaki |2002| ; iLiu &: Mirabell 120051 ). While some can be identified as supernova 



remn ants, background Active Galactic Nuclei, or faint foreground stars (e.g., I Gutierrez 
20061 ). most seem to be the result of accretion from a high-mass star onto a compact 



object. Short-term variability of some of these ULXs indicates t hat they are not simply 



unresolved superposi tions of several lower- luminosity sources (e.g.. iMatsumoto et al.l 12001 



IMBHs: 


Colbert & Mushotzkv 


Madhusudhan et al. 1 


2006h or 



super-Eddington accretion. Super-Eddington accretion can be achieved in two different 
ways: an inhomogeneous accretion-disk structure in which the photon flux is spatially sepa- 
rated from the bulk of the matter influx, or strongly anisotropic radiation. 



Anisotropic emission (IKing et al.l l200ll ) may originate, e.g., from jet sources associ- 
ated with accretion onto solar-mass black holes (microquasars) in which the jet is oriented 
at a small angle with respect to our line of sight - - the so-called "microblazar" model 



( iGeorganopoulos. Aharonian fc Kirk 1 120021 ; iKording. Falcke fc Markoff 1 120021 ) . The latter 
model is now considered unlikely bas ed on recent observations of X-ray ionizat ion of optical 
nebulae associated with some ULXs (jPakull fc Mirioni 1120031 ; iGutierrez II2006I ). 



A promising mechanism for sustaining super-Eddington accretion in a radiation-dom- 



™ 1 

Arons 


1992; 


Gammie 


1998; 


Begelman 


2001, 


2002. 


2006) 



escape along low-density regions (LDRs) aligned with predominantly vertical magnetic field 
lines. The accretion flow is concentrated in thin, optically thick high-density regions (HDRs), 
where the magnetic pressure of field lines oriented predominantly within the disk dominates 
over the gas pressure, thus confining the gas and preventing the radiation pressure from 
disrupting the accretion flow. In such a configuration, the total disk luminosity can exceed 
the Eddington luminosity by a factor approximately equal to the ratio of the magnetic field 



pressure to the gas pressure (jBegelman 



20011 . 120021 ) in the HDRs. 



Given the variety of different promising candidate models for the nature of ULXs, the 
obvious question is: Can one distinguish between these models with ULXs' X-ray spectra? 
If so, what are the characteristic spectral features of the individual models? From an ac- 
creting IMBH it is generally believed that a multi-color disk blackbody (MCDBB) spectrum 
with an inner temperature of kT in ~ 0.1 — 0.3 keV (that is, the disk temperature at the 
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innermost stable circular orbit) and possibly a high energy component from Compton up- 
scattering in a tenuous, hot corona would be a realistic phenomenological description of the 
spectru m. However, a recent, more detailed analysis of non-LTE accretion flows around 
IMBHs (jHui. Krolik & Hubenyl 120051 ) has indicated that the effects of black-hole rotation 
and Compton scattering within the disk may very well lead to much higher apparent disk 
temperatures, up to kT ~ 1 keV, in addition to deviations from conventional MCDBB 
spectra at both soft and hard X-ray energies due to metal opacity effects. 

The continuum spectrum from the photon bubble model might be dominated by the 
MCDBB spectrum emanating from the HDRs. However, this spectrum might be modi- 
fied during the radiation transport in the photon bubble cavities. Furthermore, the almost 
free-streaming radiation of the photon bubbles could be repeatedly Compton-reflected off the 
surfaces of the HDRs, potentially leading to strong fluorescence lines and/or radiative recom- 
bination edges, in addition to Compton reflection features from hard X-ray emission imping- 



ing upon the disk from radiation sources external to the disk ( IBallantyne. Turner fc Blaes 



2004 ; IBallantyne. Turner fc Young II2005I ). Such features may be observable, assuming it is 
not overwhelmed by other radiation sources, in particular the blackbody from the HDRs. 

Distinguishing scenarios by spectral modeling is currently very difficult for all but the 
highest quality ULX data sets, and spectral fit ting to these ULX spectra gives contradic- 
tory results. For example, iFeng &: Kaaretl (120051 ) recently performed a detailed spectral and 
timing analysis of archival XMM-Newton data on 28 ULXs that had sufficient photon statis- 
tics to allow for meaningful fitting with models more complicated than a simple power-law. 
They found that their continuum spectra fell into three general categories: (1) Optically thin 
bremsstrahlung-dominated, quasi-thermal spectra with temperatures of kT ~ 0.6 - 0.8 keV, 
characteristic for X-ray emission from young supernova remnants; (2) MCDBB spectra of 
temperatures kT ~ 0.1 - 0.4 keV, plus occasionally a hard X-ray power-law, as possible in 
the case of accretion onto IMBHs; (3) MCDBB spectra at temperatures kT y 1 keV, plus a 



power -law component dominating at lower energies. Based on these results, IFeng fc Kaaret 



( 120051 ) suggest that ULXs may, in fact, not be a homogen eous class of objects. This is 



in accord with recent results of iMadhusudhan et al. I (120061 ) that realistic stellar evolution 



and population synthesis calculations suggest that the expected rate of captures of mas- 
sive stars by IMBHs may not be sufficient to produce the total observed number of ULXs. 
Winter et al. I (120061 ) classified ULXs into low/hard and high/soft states based on their X- 
ray luminosities and spectra, assuming the ULXs were IMBHs accreting in states similar 
to galactic X-ray binaries. They found that the high/soft ULXs were grouped around one 
of two blackbody temperatures: one grouped around 1 keV and one grouped around 0.1 
keV, further indication that ULXs are not a homogeneous group. The long term monitoring 
of some ULXs such as NGC 5204 X-l shows that its X-ray spectrum hardens as its flux 
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increases, opposite to what has been observed in Galactic X-ray binaries ([Roberts et al 
2006h . 



Stobbart. Roberts fe Wilms I (120061 ) found that 6 out of 13 of the most-observed ULXs 



were fit approximately equally as well with a cool MCDBB and hard power-law as with a 
soft power- law and warm MCDBB. Even more surprisingly, they found that most spectra 
(10 out of 13) were fit best by a cool, kT ~ 0.2 keV blackbody and a warmer fcT, n ~ 2 keV 
MCDBB, and 11 out of 13 were fit by a cool MCDBB and an optically thick (r ~ 10) corona. 
These spectra cou ld also be explained by a hot, inner (optically thick) plasma sphere and 
a cool outer disk (lAgrawal fe Misra 1120061 ). in which case the cool blackbody temperature 
could not be directly related to the black hole mass by virtue of the relation between the 
inner disk temperature and the black hole mass. 

In summary, current diagnostics for the distinction between different models for ULXs 
indicate that ULXs may in fact not be a homogeneous set of physical objects, different models 
may apply to different sources, and in many cases, observational results are inconclusive. 
Consequently, more detailed model predictions of various ULX models might be helpful in 
identifying additional diagnostics which may be used to confirm or rule out such models. 
In this paper, we investigate detailed spectral features resulting from the photon bubble 
model, both in the continuum and possible emission line and radiative recombination edge 
features resulting from the radia tive feedback between the HDRs and LDRs. An Fe K lin e 
has been observed in M82 X-l (jStrohmayer fe Mushotzky 1 120031 ; lAgrawal fe Misra 1120 06); 
absor ption edges have been observed in the 0.1-1 keV range in M101 ULX-1 (IKont et al. 
20041 ); and possible edges at ~ .7 keV and ~ 8 keV, and an emissio n feature at ~ 6 keV has 



been seen in Holmberg IX X-l (IDewangan. Griffiths fe Rao 1 120061 ). In the ULX M51 X-26 
emission features have been seen at 1.8 keV, 3.24 keV, 4.03 keV, and 6.65 keV with Chandra 



(ITerashima fc Wilson 1 120041 ) and at 6.4 keV with XMM-Newton (IDewangan et al. 1 120051 ) . 



This indicates that the detection of emission lines and/or radiative recombination edges is 
currently feasible, at least for bright ULXs, and predictions of expected line features from 
various models will thus be useful as an additional model diagnostic. 



Studie s of X-ray reflection in inhomogeneous accretion disks have been done by lRoss. Fabian fe Ballant 



(120021 ) and iFabian et al. I (120021 ). who found that multiple Compton reflections of emission 
from a corona above an irregular or corrugated accretio n disk in Seyfer t gala xies lead to 
much stronger reflection features than single reflections. iMerloni et al. I (120061 ) considered 
Compton reflection in a radiation-dominated inhomogeneous accretion disk. Their disk was 
a low density Comptonizing plasma with high density clumps off of which radiation was 
reflected. Their simulations produced lines which may be visible in narrow line Seyfert 1 
galaxies, although detection may be difficult due to interstellar absorption. So far, inves- 
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tigations of photon bubbles in ULXs have focused on radiation-hydrodynamics (see, e.g 



Ruszkowski fe Begelman 



2003 



essent ially unexplored. iBallantyne. Turner fc Blaes 



Turner et al. 112005 ) leavi ng t he expected spectral features 



(120041 ) and lBallantyne. Turner fc Young 



( 120051 ) studyed X-ray reflection off accretion disks in the photon bubble model for Seyfert 
galaxies, but using an unspecified external irradiating hard X-ray source, such as a tenuous, 
hot corona above the accretion disk surface, while we are focusing on the self-consistent local 
radiation feedback between the HDRs and LDRs within the inhomogeneous accretion flow. 
Furthermore, those authors did not explore parameters appropriate for ULXs. 

In this paper we make predictions of the X-ray spectra from the photon bubble model 
of ULXs with a Monte Carlo/Fokker-Planck code coupled with the XSTAR program for 
X-ray reprocessing in the HDRs (reflection). We will be primarily interested in predictions 
of spectral features in addition to a warm MCDBB and soft power-law component, although 
this is not a study of any individual object. In addition we make predictions of a power- 
law component that extends above 10 keV which may not be detectable by Chandra or 
XMM-Newton but may be by Suzaku or future telescopes. 



2. Model Setup 



The simulat ion parameters were based on solutions found in iBegelman I (120011 ) and 
(120021 ). We summarize these solutions in § 12.11 and describe the simulation 



Begelman 



technique in 



2.1. Disk Distribution 



The physical picture emerging from numerical simulati ons of the photon bubble in- 
stability (jRuszkowski fc Begelman 1120031 ; iTurner et al. 1120051 ) points towards a propagating 
pattern of slab-like shock trains (the HDRs) slanted with r espect to the plane of the accretion 
disk. A solution to the disk structure has been derived by IBegelman I (120021 ). who found the 
relevant disk parameters to depend on: the accretion rate, m = Mkc/ ' AttGMbh-, ft, the ratio 
of the HDR's gas pressure to the magnetic pressure, and a, the Shakura-Sunyaev viscosity 
parameter. Here k is the opacity (assumed to be dominated by Thomson scattering; k = 
0.4 cm 2 g -1 ), c is the speed of light, G is the gravitational constant, and Mbh is the mass 
of the black hole; We used his solution to set up our simulations. Throughout this work we 
assumed /3 = 0.1 and m = M BH /M G = 10 (resulting in an Eddington luminosity of 1.3 ■ 10 39 
erg s _1 ). 
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The Eddington enhancement factor, the ratio of the disk's flux to the Eddington flux, 



is: 



/ = 



Edd 



n H PR _ 3 m£> 
" 2 5r 



n 



avg 



where n avg = ^n HDR n LDR is the geometric average of the LDR and HDR densities; r is the 
distance from the black hole in terms of gravitational radii ( R g = GM/c 2 ); S is the ratio of 



the height to the radius (5 = h/R); and D 
following equation for 5: 



(6/r) 1 ' 2 . iBegelman I (120021 ) derives the 



~ max 



0.2 



0.3 



0.1 



o.i 



4/21 



-4/13 



a 



a 



0.01 



-5/13 



(mD) 



5/13 m -l/13 r -5/26 



0.01 



-5/21 



(2) 



Note tha t in Eqn. |2 



we c orrected a typographical error in Eqn. 13 of IBegelman I (120021 ). 



Based on IBegelman I (120021 ) 's expression for the average density, and equation [TJ the following 



equations were derived for the densities in the HDR and LDR: 



n HDR = 6.0 • 10 18 (mD) 2 (ma) _1 5~V 



nLDR 



3.6 ■ 10 18 (ma)- 1 S- 2 r- 1 / 2 



-5/2 cm -3 



cm 



(3) 
(4) 



Although these densities refer to the disk midplane, the IShakura &: Sunyaevl (1l973l ) solution 
for a radiation-dominated accretion disk shows that the density will be essentially constant in 
the vertical direction. Thus, we consider these densities to be the density at the photosphere. 
The temperature for the HDRs was found by assuming that all the gravitational potential 
energy is released in the HDRs, an d thus, the flux emit t ed at the photoplane is set equal to 
the blackbody flux. This gives the IShakura fc Sunyaevl (119731 ) expression: 



kT HDR = 5.09 m" 1/4 (mD) 1/4 r- 3/4 keV. 



(5) 



The "wavelength" of the plane-parallel shocks — i.e., this distance between HDRs — is 
$iven by 

A = Pa5R. (6) 



For the proton temperature in the LDRs, the virial temperature was used: 



Ip 



GM BH m p 
R 



(7) 



where m p is the proton's mass. 
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The electron temperature in the LDR is expected to be significantly lower than the 
proton temperature due to Compton cooling by radiation from the HDRs. The LDR's 
electron temperature was determined by implicitly solving a Fokker-Planck equation with 
the MC/FP code; see § 12.21 below. The magnetic field in the LDR was determined by 
assuming it to be in equipartition with the electrons. 

The Eddington ratio is dependent on radius (equation [TJ , so for a given disk, outside 
of a certain radius {r inhom ) I < 1, and the disk becomes a homogeneous, Shakura-Sunyaev 
disk. Outside of this radius, the disk was represented by MCDBB with kT oc r -3 / 4 . 



2.2. Simulation Description 



Equations [2] through [7] were used to determine the simulation parameters at various 
radii. To find the total spectrum of the disk, simulations were run at evenly-spaced radii 
and the results were averaged, weighted by disk area. At each radius, we simulated one 
individual LDR sandwiched by two HDRs in a plane-parallel geometry (see Fig. [1]). Within 
the plane-parallel geometry the LDR is divided into 40 zones, 4 radial and 10 vertical. It 
seems r easonable to assu me, that, since A is much smaller than the radiation pressure scale 
height (IBegelman 1 120021 ) . that the fraction of photons which escape (given by the ratio of 
the escape area to the region's total surface area), will be ~ 1%. Larger escape fractions 
could lead to smaller ionization parameters and greater reflection components, however, this 
was not explored in this work. Multiple Compton reflections are certainly possible with such 
a small escape fraction, however for nearly all of our cases, the reflected component was 
completely drowned out by the blackbody; see below. As a result, further iterations did not 
produce significantly different spectra. 

The radiation transfer and self-consistent balance between electron heating and cool- 
ing within the LDR are si mulated with the two-dimensiona l Mo nte Carlo /Fokker-Planck 
(MC/FP) code described in lBottcher. Jackson fc Liang I (120031) an d Bottc her k, Liang I (1200 ll ). 



Pozdnyakov. Sobol fc Sunvaev J 19831) fo r Comp- 



This code uses the Monte Carlo method of 
ton scattering and the implicit Fokker-Planck method of iNayakshin fc Melial (119981 ) for the 
evolution of the electron distribution in a two-temperature plasma with a given proton tem- 
perature. In each zone, the Fokker-Planck equation is implicitly and independently solved. 
The Fokker-Planck technique takes into account heating/cooling by Coulomb/M0ller inter- 
actions, Compton scattering, synchrotron/cyclotron processes. 



Initially, the emission emanating from the HDRs was represented by blackbody spec- 
tra inserted at the upper and lower boundaries. Photons were also produced by syn- 
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chrotron/cyclotron processes in the LDRs. Photons were subjected to Compton scatter- 
ing and reflection off the HDRs. Escaping photons at the outer boundaries were added to 
an event file for later spectral extraction; this involves placing each escaping photon in a 
particular energy bin. 

The proton temperature in the LDR was calculated from equation [71 and the electron 
temperature was calculated numer i cally within the MC/FP simula t ion (seelNavakshin fc Melia 
19981 : IBottcher. Jackson & Liang 1 120031 : IBottcher & Liang I l200ll : iFinke k Bottcher II2005T ). 
which was run until the electron temperature reached a stable equilibrium. The MC/FP 
simulation was run to extract the photon flux and spectrum incident on the boundaries be- 
tween the HDRs and the LDR. To calculate the expected spectral features from fluorescence 
line emission, radiative recombination, and Compt on reflection, the impinging spectrum was 
used as an input into the latest version of XSTAR (jKallman fc Bautista Il200ll ). Solar abun- 
dances of the most profuse astrophysical elemen ts (H, He, C, Ca, N, O, Ne, Mg, Si, S, Ar, 
Ca and Fe) were assumed in the HDR, based on lGrevesse. Noels fc Sauval I (119961 ). XSTAR 
was run in constant pressure mode; the pressure was calculated from the specified density 
and the ideal gas law. The main parameter in determining the reflection spectrum's shape 
and intensity is the ionization parameter: 

AtxF 

t = — • (8) 

n-HDR 

Note that XSTAR uses the flux calculated between 1 and 1000 Ry (13.6 eV to 13.6 keV). In 
order to circumvent XSTAR's limitation to densities n e < 1 ■ 10 17 cm -3 , the impinging flux 
was re-scaled to keep the ionization £ at the value corresponding to the physical situation. 
The inverse flux scaling was applied to the XSTAR output spectrum. Keeping £ constant 
insures that the scaling of a flux dominated by recombination features, F oc n e is properly 
recovered. Testing of XSTAR with various densities but the same £ seems shows that this 
is appropriate for the energy range of interest. The resulting reflection spectrum was then 
added to the intrinsic blackbody from the HDR as boundary sources in a second run of the 
MC/FP code for the final evaluation of the emanating X-ray spectrum. For an example of 
a simulation, with spectra at different radii, see Fig. [2j 



3. Results 

The simulation parameters and fit results can be seen in Table [TJ The important disk 
parameters and are plotted as a function of radius in Figs. [3] to 

The transition between the expressions for 5 (Eqn. [2]) can be readily seen in panel (a) of 
Fig. [3X with the first expression being used at lower radii. At higher a, the transition occurs 
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at lower radii, or not at all. Note that several simulations violate the thin-disk condition 
(5 < 1), especially the high rh and/or low a ones. These solutions should thus be viewed 
with caution. 

The ionization parameter decreases with radius, as is observed in our simulations, seen 
in panel (b) of Figs. [3] to [51 At no point does £ drop below ~ 2000 erg cm s _1 , so the 
reflection component will be small (< 12% of the total) for all cases. 

The optical depth, r, is proportional to Xtildrj and is essentially always increasing with 
r. Thus the spectra have more luminosity from Comptonization at larger radii. It is also 
larger for smaller accretion rates due to its 6 dependence (r oc Xuldr oc Suldr oc <5 _1 ) 
by virtue of 5 increasing with increasing rh. Thus we expect a greater Comptonization 
component for smaller rh. We can see that for our simulations, this is in fact the case (Fig. 
E]). We can also see in Tabled] that the fractional power emitted through Comptonization 
decreases with increasing rh, as expected. This also corresponds to the decrease in T, the 
photon index (ELe oc E~ t+2 ) with increasing rh. One can also see in Figs. O— [5] that as a 
increases, r increases. Again, this is because r oc <5 -1 , and as a increases, 5 decreases. 

The ratio I peaks at r ~ 12, as can be seen in panel (d) of Figs. [3]-[5j This is close to 
where the peak flux is expected in a typical Shakura-Sunyaev accretion disk. Also note that 
I increases with increasing a, as one would expect from eqn. [1] (I oc S^ 1 ). 

For each simulation, the total reflection component from XSTAR was summed up for 
all radii. The soft X-ray spectra were fit with a MCDBB. The MCDBB and reflection com- 
ponents were subtracted from the total spectra revealing the Comptonization components. 
This component was fit with a power-law and exponential cutoff (ELe oc E~ t+2 e~ E ' Ecut °ff) 
above 25 keV. The results of this decomposition and fits are summarized in Table (TJ the in- 
dividual fit components for the a = 0.01 simulations can be seen in Fig. [3 Note that these 
fit components are the sum of the simulations at all radii. The only simulation that has a 
significant reflection feature visible in its spectrum is Sim. 2, with a blended Fe XXV/XXVI 
feature at ~ 9 keV with an equivalent width of ~ 640 eV, and an edge at ~ 0.87 keV (O K). 
This simulation also has the highest fraction of its emission from the reflection component. 
One can see in Fig. [2] that the reflection component is strongest at r = 120 — 280. 

In Fig. [2j one can see that the Comptonization component gets stronger at larger radii, 
as one would expect with the increasing optical depth. The overall spectra are generally 
dominated by the high-energy Comptonization component, especially for higher rh. The 
Comptonization component in many cases extends far down in photon energy into the ~ 
few keV range and often dominates the total bolometric luminosity, even when the total 
spectrum appears to be dominated by a thermal soft X-ray component. 
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This component contributes to the spectra even at low energies, down to ~ 1 keV, as seen 
in Fig. [7J In Fig. [2J one can see that the Comptonization component gets stronger at larger 
radii, as one would expect with the increasing optical depth. This explains the significant 
Comptonization component to the lower energy part of the spectrum. Eventually, though, 
one gets to the radius where photon bubbles can no longer be sustained (r in hom), and there 
is no longer a Comptonization component, as seen in r > 520 in Fig. [21 Results for other 
simulations are similar to Fig. [2j As a increases, the r inhom increases, as can be seen in Figs. 
EH51 and Table □ 

Spectral pivoting can be seen in the all of the simulations in Fig. [61 This is due to 
the approximately same amount of luminosity in the Compton components of the different 
simulations. The energy dissipated in the LDRs is related to the proton temperature, which 
in our simulations does not depend on accretion rate. This is not an unrealistic assumption, 
as one would not expect the protons to radiate significantly. However, when the LDRs 
become more optically thick the energy dissipation becomes more efficient due to multiple 
Compton scatterings. 



4. Discussion 

Our simulations explore 0.5-10 keV luminosities up to ~ 8 • 10 39 erg s _1 which does not 
include some of the brightest ULXs. For our setup, luminosities much greater than this will 
violate the thin disk condition, except possibly for the highest value of a. It does not seem 
likely that photon bubbles could support much greater luminosities. 



Our results are significantly different from previous, similar works. iMerloni et al. 



( 120061 ) perform similar calculations for a low density accretion region with imbedded op- 
tically thick clumps and considered reflection between these clumps. They used ionization 
parameters comparable to ours (~ 3000 erg cm s _1 , among others) and achieved much 
stronger reflection components with no ~ 9 keV feature. This is a consequence of the sub- 
stantially lower temperature of the soft blackbody (ultraviolet) photons used as an input 
into the Comptonization scheme in their model: their Comptonization component extends 
down to < 0.1 keV energies (in contrast to ~ a few keV in our case). For comparable £, 
calculated from the 1 — 1000 Ry flux the number of photons above the Fe K edge is thus much 
smaller than in our case. Therefore, heavy metals are almost fully ionized in most of our 
simulations, leading to weak or absent line featu res and a relatively weak Compton reflection 



component. iBallantyne. Turner &: Blaes I (12004 ) explored reflections from a hot corona off of 



an inhomogeneous disk and found the spectrum can differ significantly from reflection off a 
homogeneous disk. They also had stronger reflection features than our results, due to lower 
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ionization parameters, and no ~ 9 keV feature. 



Feng fc Kaaretl fl2005f ) and lStobbart. Roberts fc Wilms I fl2006f ) fit ULXs' X-ray spectra 



with various models; when fit with a MCDBB or MCDBB and low-energy blackbody, their 
fits had inner disk temperatures similar to ours. However, we find that our model is unable 
to explain a soft-excess. We note that no ~ 9 keV feature has b een detected in the spectra 



of ULXs, although features at other energies have been detected ( Stroh maver fc Mushotzky . 
20031 : lAgrawal fc Misra II2OO6I : iKont et al. hooi bewangan. Griffiths fc Rao Il2006h . although 



these features could originate from a wind rather than reflection features in the disk. ULX 
observations above ~ 10 keV by, Suzaku or the next generation of hard X-ray imaging in- 
struments might be able to detect the hard power-laws we predict, although we realize such 
observations would be difficult. Suzaku observations of two ULXs in NGC 1313 did not 



detect any component above 10 keV above the background (IMizuno et al. 1 120071 ). although 
their MCDBB + power-law fits do seem to agree with our spectra, and the variability of 
these sources is much less than observed in Galactic black hole candidates. 

Unfortunately, our simulations do not reproduce a soft excess as has been observed in 
many ULXs. It is possi ble that the soft excess originates from Compton downscattering of 
radiation by wind (e.g., iBegelman Il200ll ). and so a lack of its production in our simulations 



does not disprove the photon bubble model for ULXs. 

If the excess is explained, our model could explain ULXs well fit with a ~ 1 keV 
blackbody. Higher energy observations (> 10 keV) could determine this; if they found a 
stable hard X-ray power-law V ~ 2 — 2.4 for a long period of time, this would be evidence 
for the photon bubble mo del in ULXs, due to the fact that a lower mass compact object 
would have less variability (IKalogera et al.ll2004l ) . A ~ 9 ke V feature could also be considered 
a signature of photon bubbles in ULXs; it has not been found in any other simulation of 
accreting black holes. 

We no te that our spectra are similar to the ve ry high state of X- ray binaries such as 
GX 339-4 jBelloni et al. Il2006h and GRO J 1655-40 JSaito et al. Ibood ). which have similar 
photon indices as found in our simulations. For smaller accretion rates, the photon bubble 
model may be a viable model to explain the very high state of X-ray binaries. 

Since this is an exploratory study, we have made a number of simplifying assump- 
tions. We have assumed a plane-parallel geometry; zones could be corrugated or vary in 
random ways, which could lead to multiple reflections, increasing the reflection component 
(IFabian et al. Il2002l ). We have assumed constant densities across LDRs. This may alter the 
Comptonizing region of the spectra, but probably would not lead to significant differences. 
We have completely neglected several items. Taking into account General Relativistic effects 
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could gravitationally broaden the spectrum. However, most of the emission originates at 
r ~ 12, too far out to be greatly affected by GR. We have assumed n HDR and n LDR to be 
constant in the vertical direction, neglecting; c olor c orrections; this may be validated by the 
more detailed vertical solution of iBegelman I (120061 ). who found that the corrections should 
be minor. Time variability is beyond the scope of this study and is poorly understood, 
and thus would involve poorly constrained parameters. It is also possible that the spectra 
could be modified by further Comptonization in a corona above the disk, which we have 
also neglected. This would make the spectra even harder than they already are, possibly 
too hard to match observations. Ho wever, recent magnetohydrody namic simulations have 
shown difficulties in creating coronae (IHirose. Krolik. fe Stone 1120061 ) and the photon bubble 
model may be an alternative to the standard disk-coron a geometry. Futu re simulations could 



take advantage of more detailed analy tic solutions (e.g.. IBegelman 1120061 ) and hydrodynamic 
simulations (e.g. 
disk wind. 



Turner et al. 1120051 ) and could include reprocessing in a corona and/or a 
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Table 1. Simulation and fit parameters, m is the accretion rate, a is the Shakura-Sunyaev viscosity parameter, kT in 
is the inner disk temperature, and r inhom is the radius above which photon bubbles cannot exist (I < 1). Results of 

our simulations are: total luminosity over the Eddington luminosity is given (L tot /LEdd)\ the 0.2 — 10 keV luminosity 
over the Eddington luminosity (Lx/ LEdd); the fitted photon index, T; the high energy exponential cutoff; and the 

fraction of the total luminosity in the MCDBB component (Jbb), m the Comptonization component (fcompt), and the 
reflection component {frefl)- Note that these do not all add up to unity due to rounding. 



Sim. No. 


rh 


a 


kT m [keV] 


^inhom 


Ltot/L Edd 


Lx/LEdd 


r 


Ecutoff [keV] 


/ss 


fCompt 


frefl 


1 


200 


0.01 


1.4 


280 


2.8 


1.3 


2.18 


460 


0.22 


0.73 


0.05 


2 


400 


0.01 


1.6 


472 


4.4 


2.5 


2.42 


440 


0.26 


0.63 


0.11 


3 


900 


0.01 


2.0 


860 


6.9 


4.7 


2.42 


290 


0.56 


0.40 


0.05 


4 


200 


0.1 


1.4 


482 


5.0 


1.7 


1.96 


280 


0.17 


0.82 


0.01 


5 


400 


0.1 


1.6 


804 


6.5 


2.9 


2.12 


220 


0.27 


0.74 


0.01 


6 


900 


0.1 


2.0 


1456 


9.4 


5.3 


2.24 


140 


0.39 


0.59 


0.02 


7 


200 


0.5 


1.4 


698 


5.8 


1.8 


1.75 


190 


0.14 


0.85 


0.01 


8 


400 


0.5 


1.6 


1164 


8.5 


3.1 


1.78 


160 


0.17 


0.79 


0.05 


9 


900 


0.5 


2.0 


2100 


11.9 


5.6 


1.95 


160 


0.31 


0.68 


0.02 
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Fig. 1. — The simulation geometry. Each simulation is of one part of the disk. The LDR 
is divided into 40 zones. Afterwords, simulations of different parts of the disk are averaged, 
weighted by area. See the electronic edition of the Journal for a color version of this figure. 
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Fig. 2. — Spectra at various radii and the total spectrum for Simulation 2 
edition of the Journal for a color version of this figure. 



See the electronic 




Fig. 3. — The parameters S (a), £ (b), r (c), and I (d) as a function of radius, for m = 200 
(solid black), m = 400 (dotted red), and m = 900 (dashed green). These simulations have 
a = 0.01. See the electronic edition of the Journal for a color version of this figure. 
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Fig. 4. — Same as Fig. [3] except for a = 0.1. See the electronic edition of the Journal for a 
color version of this figure. 
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Fig. 5. — Same as Fig. [3] except for a = 0.5. See the electronic edition of the Journal for a 
color version of this figure. 
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Fig. 6. — Total spectra for (a) a = 0.01, (b) a = 0.1, and (c) a = 0.5. For all graphs we 
plot m = 200 (solid black), m = 400 (dotted red), and m = 900 (dashed green). See the 
electronic edition of the Journal for a color version of this figure. 
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Fig. 7. — The spectra broken into components for a = 0.01 and (a) rh = 200, (b) rh = 400, 
and (c) rh = 900. The solid black line is the reflected component, the dotted red line is the 
MCDBB, the dashed green line is the Comptonization component, and the long dashed blue 
line is the total spectrum. See the electronic edition of the Journal for a color version of this 
figure. 



